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Abstract 



We study the theoretical accuracy of various methods that have been proposed to 
measure the luminosity of the LHC pp collider, as well as for Run II of the Tevatron pp 
collider. In particular we consider methods based on (i) the total and forward elastic data, 
(ii) lepton-pair production and (iii) W and Z production. 



1 Introduction 



The Large Hadron Collider (LHC) being constructed at CERN, will generate proton-proton 
collisions with total cm. energy of 14 TeV with a design luminosity L = 10^^ cm~^ s~^. 
The experiments at this new facility will have a high potential to discover New Physics and 
to make various precision measurements, see, for example ||. Both the general purpose 
pp experiments, ATLAS and CMS |Q, will provide high statistics data samples, and the 
accuracy of the precision measurements will be limited by systematic effects and, in many 
cases, by the uncertainty in the measurement of the luminosity C. For example, precision 
measurements in the Higgs sector of typical accuracy of about 7%, for a wide range of possible 
Higgs mass, require the uncertainty in the luminosity to be < 5% 0. Both LHC experiments 
are aiming at measuring the luminosity to 5%. Recall also that the enhancement in luminosity 
that will be achieved in Run II of the Fermilab Tevatron pp - collider will herald a new era of 
precision studies P, 0- 

An obvious requirement for the success of these precision measurements is that the uncer- 
tainties in the theoretical calculations of the cross sections for the basic processes, which are to 
be used to determine the luminosity, should match the desired experimental accuracies. 

In general, there are two possibilities to determine the luminosity - either (i) to measure a 
pair of cross sections which are connected quadratically with each other, or (ii) to measure a 
cross section whose value is well known or which may be calculated with good accuracy. The 
well-known example of the first possibility is the measurement of the total Utot and differential 
forward elastic cross sections which are related by the optical theorem; see, for example, |^ 
for recent experimental discussions. This method is discussed further in Section 2. 

Two types of processes stand out as examples of the second possibility to measure the 
luminosity. First there is exclusive lepton-pair production via photon-photon fusion 

7+7- n^ 
P P ^ P + ri + P (1) 

where I = e or fi. To the best of our knowledge this proposal originated in 0. A luminometer 
for the LHC based on measuring forward e^e~ pairs (of invariant mass Mee < 10 — 20 MeV and 
transverse momentum pt{ee) < 10 — 20 MeV) was proposed in |T^, while Ref. |[ll|] concerns 
the central production of fi~^fi' pairs (with (M^^) ~ 20 GeV and PtijJ'lj) ~ 10-50 MeV). 
Lepton-pair production is the subject of Section 3. 

Nowadays attention has also focussed on W and Z production as a possible luminosity 
monitor, both for Run II at the Tevatron and for the LHC see, for example, [|T^. The reason is 
that the signal is clean, and the production cross sections are large and can now be calculated 
with considerable theoretical accuracy. We discuss this possibility further in Section 4. 

In principle, we may monitor the luminosity using any process, with a significant cross 
section, which is straightforward to detect. For example, it could be single- (or two-) pion 
inclusive production in some rapidity and pt domain (for example = 5 — 10 GeV; |?7| < 1) or 
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inclusive n~^jj,^ production in a well-defined kinematic domain, etc. In this way we may control 
the relative luminosity and then calibrate the "monitor" by comparing the number of events 
detected for the "monitor" reaction with the number of events observed for a process with a 
cross section which is already known, or which may be calculated with sufficient accuracy. 

For some applications better accuracy can be achieved by measuring the parton-parton 
luminosity. A discussion is given in Section 5. Finally, Section 6 contains our conclusions. 



2 Elastic scattering 

First we discuss the classic method to measure the luminosity, that is using the observed total 
and forward elastic pp (or pp) event rates. Neglecting Coulomb effect^, the elastic cross section 
is given by 



da, 



el 



dt 



2 

tot /I I 2 



! + /)• (2) 



The ratio p of the real to the imaginary parts of the forward amplitude is small at Tevatron-LHC 
energies and can be estimated via dispersion relation techniques 

_ TT glnatot(g) . . 

^ - 2 dins • 

For example, at the LHC energy, y/s = 14 TeV, it is predicted to be p = 0.10 — 0.12; for a 
recent estimate, see, for example, Ref. |T^. Even the largest conceivable uncertainty in p, that 
is Ap = ±0.02, leads only to an uncertainty of less than 0.5% in daei/dt. Thus if we measure 
the number of events corresponding to the elastic scattering and to the total cross section, then 
we may determine both the luminosity C and cxtot (since iVei oc a^C, whereas A^tot oc <JtC). The 
main problem in that it is extremely difficult to make measurements at the LHC in the very 
forward region, and so it is necessary to extrapolate elastic data from, say, \t\ > 0.01 GeV^ to 
t = 0. 

On the other hand it was found at the ISR that the observed 'local' slope 

Bit) ^ dHd^ 



depends on t. In particular at ^/s = 62 GeV ||T6 



AB = B{0) - B{\t\ = 0.2 GeV) ~ 2 GeV"'. (5) 

^Strictly speaking, we have to use Coulomb wave functions, rather than plane waves, for the in- and out- 
states for elastic scattering between electrically charged protons. To account for this effect, the strong interaction 
amplitude should be multiplied by the well-known Bethe phase ||l3|] {Ad — > exp{iSB)Aci). The uncertainty 
in Bethe phase Sb due to the non-point-like structure (electric charge distribution) of the proton is about 
a — 1/137. It leads to an uncertainty of about ±0.01 in the Re/Im ratio, and hence to less than a 0.1% 
correction for the imaginary part of the strong amplitude ^oi- A recent discussion of the Coulomb phase for 
t ^ can be found in ||l3| . 
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Note that if this t dependence were to be solely due to a non-linear Pomeron trajectory | T7[ 



then the difference would increase as the logarithm of the energy (Ins); giving, for example, 
AB ~ 4.8 GeV~^ at the LHC energy. The t-dependence of the local slope B(t) has recently 
been determined using a model which incorporates all the main features of high energy soft 
diffraction. That is the model embodies 



(a) the pion-loop contribution to the Pomeron pole (which is the main source of the non- 
linearity of the Pomeron trajectory ajp(t)), 

(b) a two-channel eikonal to include the Pomeron cuts which are generated by elastic and 
quasi-elastic (with A^* intermediate states) s-channel unitarity, 

(c) the effects of high-mass diffractive dissociation. 



The parameters of the model are a(0) and a' of the bare Pomeron trajectory, two param- 
eters describing the elastic proton-Pomeron vertex, and the triple-Pomeron coupling and its 
slope |]T5| . The values of the parameters were tuned to describe the observed pp (or pp) elas- 
tic differential cross sections throughout the ISR-Tevatron range, restricting the description 
to the forward region \t\ < 0.5 GeV^. Note that there are two effects, leading to a t depen- 
dence of B, which act in opposite directions. First, the non-linear pion-loop contributions 
to the Pomeron trajectory, Aa;jp(t), lead to a contribution AB]p{t) = Aa;jp(t)lns in single 
Pomeron exchange, which increases as |t| 0. On the other hand the absorptive (rescattering) 
corrections, associated with eikonalization, lead to a dip in daei/dt (whose position moves to 
smaller |t| as the collider energy, y/s, increases), with the result that the local slope B{t) grows 
as —t approaches the position of the diffractive minimum; that is B{t) decreases as \t\ — > 0. 
Fortunately, just at the LHC energy, these two effects almost compensate each other. 

So far it has proved impossible to sum up in a consistent way all the multi-Pomeron dia- 
grams. For this reason two versions of the model were studied in Ref. |jl5| with maximal and 
minimal contributions from high-mass diffractive dissociation. After tuning the values of the 
parameters to describe the atot and the forward d(7e\/dt data, the two versions of the model 
predict different values of the slope at t = 0,5(0). The maximal choice gives a larger value 
B(0) = 21.9 GeV^^, while for the minimal case we obtain 5(0) = 20.3 GeV^l However in both 
cases the variation of B{t) satisfies AB < 0.2 GeV"^ for \t\ < 0.1 GeV^ or AB < 0.1 GeV"^ 
for the more restricted interval < |t| < 0.02 GeVl It means that we may neglect such a 
small variation of B{t) with t and safely extrapolate the observed cross section da^x/dt down 
to t = 0, using the slope measured, say, in the interval 0.01 < \t\ < 0.05 GeV^ The error, due 
to the variation of B with t, is expected to be less than AB ■ \t\ < 0.2%. 

In the very forward region, |t| ~ 0.01 - 0.02 GeV^ it may only be possible to measure 
the elastic cross section in an LHC run with a rather low luminosity. For the high luminosity 
runs one would then have to use another monitor. One possibility is to choose a single particle 
inclusive process, with a significant cross section, and "calibrate" it in the low luminosity run 
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by comparing with the total and elastic cross sections discussed above. It would be even better 
if we were able to calibrate several different reactions. Then to use, as a luminosity monitor, 
the reaction with the closest topology (or kinematic configuration) to the process that we are 
to study in the high luminosity run. 



3 Lepton pair production as a luminometer 

At first sight any QED process, with sufficient event rate, may be used as a luminosity monitor. 
Up to ISR energies the luminosity was measured via Coulomb elastic pp scattering, where at 
very small \t\ the influence of the dominant single photon exchange contribution is evident. Un- 
fortunately, in the LHC environment, the Coulomb interference region, |t| < (1 — 3) x 10^'^ GeV^ 
is not likely to be experimentally accessible. 

Other possibilities to consider are A isobar Coulomb excitation, {pp A"'"p), or vr^ and 
f] production, which are all mediated by photon exchange, as shown in Figs. 1(a) and 2(a) 
respectively. Surprisingly, the electromagnetic widths A"*" —>■ py, tt*^ —>■ 77 and r/ ^ 77 are 
known only to 7-10% accuracy [113]. Moreover, strong interaction effects in the initial and final 
states give non- negligible corrections, see Figs. 1(b) and 2(b). 

Lepton pair production looks much more promising as a luminosity monitor. The Born 
amplitude of Fig. 3(a) may be calculated within pure QED (see, for example, the reviews in [p!9| ), 
and there are no strong interactions involving the leptons in the final state. The only question is 
the size of the absorptive corrections arising from inelastic proton-proton rescattering, sketched 
in Fig. 3(b). Fortunately the rescattering correction is suppressed for two reasons. First, the 
main part of the Born cross section (Fig. 3(a)) comes from the peripheral region with large 
impact parameter bt, where the strong amplitude A{s, bt) is small. Second, even in the small bt 
domain, the rescattering correction is greatly suppressed due to the angular integration. 

Before we discuss the last point, we note that, in practice, it is difficult to exclude contri- 
butions coming from the reactions 

pp ^ X + l^l^ + p 

pp^X + l+l- + Y (6) 

where X and Y are baryon excitations, that is A^* or A isobars. Of course the matrix ele- 
ments of the corresponding processes (such as Fig. 3(c)) can, in principle, be determined from 
photoproduction and deep inelastic data, and may be taken into account. However the matrix 
elements of Fig. 3(c) are not known to sufficient accuracy and it is better to suppress the con- 
tributions of (PD by experimental cuts. The procedure is as follows. Recall that, due to gauge 



invariance, inelastic vertices of the type p — X -|- 7 vanish like |]T9 



V{p^X^) oc qu (7) 
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as the photon transverse momentum qu ^ 0. Unfortunately it is difficult to measure a leading 
proton with very small transverse momentum, that is qt ^ 10 — 30 MeV. So to take advantage of 
the behaviour of it was proposed |TI[] to select events with very small transverse momentum 
of the lepton pair 

Pt = I qi+t + Qi-t I < 10 - 30MeV, (8) 

in order to suppress A^* and A production. In such a case the integral over the transverse 
momenta of the photons in the Born cross section of Fig. 3(a) takes the form 

where the dominant contributions come from the regions qu ^ Pt- The values of the longitudinal 
components are 

— \ti,min\ •^i'^py (10) 

where Xi are the fractions of the momenta of the incoming protons carried by the l^l^ pair, 
and nip is the mass of the proton. 



3.1 fi^fi production 

To identify muons (and to separate them from vr^ mesons) it is necessary to observe charged 
particles after they have transversed a thick (iron) absorber. It means that we consider only 
muons which have a rather large transverse energy Et > 5 GeV. However, it is still possible to 
select events where the sum of their transverse momenta is small, pt < 30 MeV. In this case 



(2 - 3)Et/yG ~ 10-^ U ^ 10-^GeV^ (11) 



at the LHC energy. In such a configuration the main contribution to the integral @ comes 
from the domains qu ^ Pt, ^2 ^ ^ Pt q2t ^ Pt, Qu Pt- Performing the qf^ 

integrations in these two domains gives \n{pf/t2) and ln{pf/ti) respectively, and so the Born 
cross section behaves as 

da 1 ( pf ^ 
oc ^ In — 

dpi pi \t1t2 



2 « -2 IH fr) (12) 



In addition to the use of the small pt cut in order to separate the elastic process (|l|) from 
the excitation process (P), we can also exploit the different kinematics of the processes. For 



example in [11 1 it was proposed to fit the observed distribution in the muon acoplanarity angle 



(j) in order to isolate the elastic mechanism via its prominent peak at = 0. 
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3.2 Absorptive corrections to production 

To determine the effect of the absorptive, or re-scattering, correction, we calculate diagram 3(b) 
with an extra loop integration over the momentum Q transferred via the strong interaction 
amplitude shown by the "blob". The relative size of the correction, to the amplitude for 
lepton-pair production of Fig. 3(a), is 

^ _ {(Jinei/8n'^) I d'^Qtd'^qitd'^Ptcr{qi,q2;qi + Q,q2 + Q) F{Q) 
J d^qud'^Pta{qi,q2;qi,q2) F{0) 

where F{Q) represents the collective effect of the appropriate proton electromagnetic form 
factors Fn and the dependence of the strong interaction amplitude 

F{Q) = F^{ql) FM) FnHQ + qif) F^{{Q + g^)') e'Q'/\ (14) 

h is the slope of the pp elastic cross section, that is dae\/dt oc exp(6t). We can interpret the 
correction S in terms of 5* = 1 — 5, where S is the "survival probability amplitude" that the 
secondary hadrons produced in the soft rescattering do not accompany £~^£~ production. Of 
course this effect is due to the inelastic strong interaction only. Consider the hypothetical case 
with pure elastic pp rescattering (at fixed impact parameter bt, that is for fixed partial wave 
£ = bt\^/2). Then elastic rescattering only changes the phase of the QED matrix element, 
A4 Ai exp{2i(j)i), and does not alter the QED cross section a = |A^p. Thus we have 
to account only for the inelasticity of the strong interaction, that is 5 oc ai^ei as in ([T3|) p|. 
Throughout this paper we do not discuss the effects of "pile-up" events. Thus there will be 
an additional factor [W < 1) which represents the probability not to include events with 
extra secondaries coming from an almost simultaneous interaction of another pair of protons. 
Experimentally the depletion due to such pile-up events may be overcome, for example, by 
cleanly observing the vertex of fi~^fi~ production. 

The cross section a in ([13|) plays the role of the cross section for the QED subprocess of 
77 scattering through the lepton box. For the absorptive process, the photons have momenta 
qi + Q and q2 + Q in the "left" amplitude A shown in Fig. 3(b), and qi,q2 in the "right" 
amplitude A*, which is not shown. Here we have assigned the absorptive effect to A. Actually 
the numerator of (|T^) , and the symbolic^ diagram 3(b), summarize a set of Feynman graphs 
which describe the effects of the strong interaction between the protons in the amplitude A. 
Then we have to add the equivalent set of diagrams in which the rescattering effects occur in the 

^ An alternative, and more formal, explanation of 5 oc a-mei is that we have to sum up multiple pp rescattering. 
On resumming all the eikonal graphs, it turns out that one obtains ([l^) with crtot — <^ci = Cinei- 

''indeed Fig. 3(b) is truly symbolic and must not be viewed literally. The problem is that the strong 
interaction is not mediated by a point-like object, but must rather be viewed as a multipheral or gluon ladder 
(Pomeron) exchanged between the protons. It turns out that the dominant contribution comes from four 
different configurations, corresponding to the photons being emitted either before or inside the Pomeron ladder. 
Due to the conservation of the electromagnetic current, the sum of all four contributions is embodied in ( p^ 
— the only qualification is that when the photon is emitted inside the ladder, the form factor may have 
different behaviour at large Qt- In our case, when Qt is very small, this effect is negligible. 



6 



amplitude A*. The total correction to (|T^) is therefore 2Sda/dpf. The rescattering correction 
is due to the imaginary part of the strong amplitude only. The effects coming from the real 
part in A cancel with those in A*. 

At first sight the largest absorptive effects, (|13|), appear to come from the largest values of 
Qt < where Rp is the proton radius; higher values of Qt are cut-off by the proton form 

factor. However when Qt ^ qu, q2t an interesting suppression occurs. In this domain 77 
occurs dominantly in the = two-photon state. The projection of the orbital angular 
momentum on the longitudinal {z) axis is clearly zero, while the spin (which originates from 
the photon polarisations) is described by the tensor Q^tQut- After the azimuthal integration, 
the tensor takes the form leading to = 0. If we neglect the mass of the lepton, 

nil ^ Et, and higher order QED corrections, then the Born amplitude, with = 0, vanishe^. 
This is a well known result, see, for example, [Q. As shown in Ref. |2^, the physical origin of 



the suppression is related to the symmetry properties of the Jz = Born amplitude. 

The "cross section" a{qi, ^2; Qi+Q, (I2+Q) for our QED subprocess depends on the transverse 
momenta of four virtual photons. The full expression for a is rather complicated. However in 
the limit Et ^ (lu,<l2t,Qt,nT'i, we can integrate over the direction of the lepton transverse 
momentum kt (with |fct| ^ Et), and the formula reduces to the simple form 

cosh(A?7) [{Q + q2,qi)iQ + qi,q2) + (Q + qi,qi)iQ + q2,q2) - iqi,q2)iQ + quQ + q2)] 



a oc 



Efcosh'ilAr^) qhUQ + q.fiQ + q2f 

(15) 

where the notation (fci, ^2) is used for the scalar product of two transverse vectors, and where 
the rapidity difference 

A7]=\7]{t) -7]{l-)\. (16) 

It is interesting to note that the dependence of the cross section a on git, g2t and Qt follows 
directly from the symmetry properties of the process. First, gauge invariance implies that a 
must vanish when the transverse momentum of any photon goes to zero. That is a must contain 
a factor qitq2t{Q + qi)t{Q + q2)t- Next, it has to be symmetric under the interchanges qi ^ q2 
and qi (Q + qi)- Finally, as discussed above, in the limit Qt ^ qu, q2t, it must vanish after the 
azimuthal angular integration of Qt has been performed. The only possibility to satisfy these 
conditions, to lowest order in Qt/Et and qu/ Et, is given by the expression in square brackets in 
(p!5|). To obtain the Arj and Et behaviour, shown in the first factor in a of ([T5|), it is sufficient 
to put Qt = and to recall the well known pure QED cross section^. 

Recall that for fi^fx^ production we select events with small transverse momentum of the 
lepton pair, pt. Now if Pt = Qit — q2t — 0, that is if Qt ^ Pt, then the last factor in expression 
(|T3|) of the absorptive cross section simplifies to 

^ ^ m + quqir-qfiQ + qi?] 

cr oc 7j , (17) 

[qKQ + qi)f 



^An analogous cancellation, which occurs on integration over the azimuthal angle, was observed in QCD for 
light quark-pair production pp ^ p + qq + p hy Pumplin see also pl| 

^We thank A.G. Shuvaev for using REDUCE to explicitly check form (nS) for a. 
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which, upon integration over the azimuthal angle of the vector {Q + qi)t, gives zero. Therefore 
the main contribution to (0) comes from values of Qt ^ Pt, and the absorptive correction 



5 - ^ C. (18) 

The coefficient C is numerically small for two reasons. First, the integration over the loop 
momentum Qt kills the logarithmic factor, ln(p^/tit2), which enhanced the original cross section 
([T2|) in the absence of absorptive corrections. To see this, note that the log came from the 
Qidqt/qf integrations, for both qu and q2t, of (|^), see (p^. However when Qt ^ the integrands 
are of the form ■ [Q^ + qt)/qt\Qt + there is no logarithmic singularity for either 

— > or IQ^ + q^l — > 0. Second, there is a suppression of C from cancellations which occur 
after integration over the remaining azimuthal angles. For example, if the invariant mass of 
the /i^/i^ pair produced at zero rapidity is M^^ = 20 GeV, then the value of the coefficient 
C = 0.14,0.13,0.09 and 0.08 for pt = 5,10,30 and 50 MeV respectively. This leads to a 
negligible correction to the cross section ([T2|) , for example 

SOmb 

25 < —— pI C < 0.02% (0.13%) (19) 

for Pt = 10 MeV (30 MeV). 

Note that the rescattering contribution of Fig. 3(b) is less singular as qu, q2t than the 
pure QED term of Fig. 3(a). Therefore the rescattering correction does not induce a sharp 
peak at = in the muon acoplanarity distribution. Thus removing the excitation processes 
(H) by fitting the distribution, will automatically suppress the rescattering correction. 



3.3 e+e production 

Unlike fi~^fi~ production, for e~^e~ production we do not need to select events with electrons of 
large transverse momentum pet- Providing the energy of the electrons ~ 5 GeV, we may use, 
say, an electromagnetic calorimeter to identify them. Therefore we may consider the small Pet 
domain where the e"'"e~ production cross section is much larger. The dominant contribution to 
the cross section comes from the region of pet of the order of the electron mass {pet ^ 1 MeV). 
As Pet is small, we cannot neglect Qt in the t channel electron propogator. As a consequence 
the integral over Qt becomes super- convergent 



d'Qt 


Qt + Qit\ 


\Qt + Q2t\ 




{Q + qi? {Q + q2)H\Q 


+ qi+Pe 


P + ml) 



(20) 



It is informative to note the origin of this form. Due to gauge invariance, the numerator vanishes 
when the photon transverse momenta \Qt + qit\ ~^ 0- The first two factors in the denominator 
arise from the photon propagators. Finally we have the factor due to the electron propagator, 
with Q ~ (0;Qj,0). Contrary to production, where the lepton propagator was driven 

by the large pt — Et of the muon, here the momentum Qt is not negigible in comparison with 
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Pet and the dominant contribution comes from the region ~ < 10 ^ GeV^. Hence the 
expected absorptive correction due to strong interaction rescattering is 

26^^ Q\ ~ ^ lO-^GeV^ < IQ-^ (21) 
47r 47r 



3.4 Conclusion on l^l production as a luminosity monitor 

We see that both and e+e~ pair production may be used as monitor processes in high 

luminosity LHC colhsions. However we must employ special selection cuts on the data. 
To suppress A^* contamination we impose small transverse momentum pt of the lepton pair, 
moreover for e^e~ production we require the individual electron to be small. We also require 
the leptons to be sufficiently energetic in order to identify them. Then the cross section for 
production at the LHC (pp — > can be calculated within pure QED and, importantly, 

we may neglect the strong rescattering effects between the protons up to 10^^ accuracy. Of 
course by requiring the muons to have high Et we have a rather small cross section. However 
here we may trace back the muon tracks and determine the interaction vertex, and hence isolate 
the interaction in pile-up events. So, in principle, fi~^fi~ pair production, with high muons, 
may be used as a luminometer in very high luminosity LHC runs. 



4 W and Z production as a luminosity monitor 

W and Z production in high energy pp and pp collisions have clean signatures through their 
leptonic decay modes, W —>■ lu and Z —>■ , and so may be considered as potential luminosity 
monitors |T2|. A vital ingredient is the accuracy to which the cross sections for W and Z 
production can be theoretically calculated. The cross sections depend on parton distributions, 
especially the quark densities, in a kinematic region where they are believed to be reliably 



known. Recent determinations of the W and Z cross sections can be found in Refs. ^5|, |26 
Here the situation has improved, and next-to-next-to-leading order (NNLO) predictions have 
been made|. The most up-to-date values are reproduced in Fig. 4 and Figs. 5,6 [^, 

To estimate the accuracy with which the W and Z cross sections are known, we start with 
Fig. 4. Fig. 4 was obtained using parton distributions found in LO, NLO and NNLO global 



analyses of the same data set pH]. It is relevant to summarize the contents of these plots. The 



predictions labelled LO, NLO and NNLO can be schematically written as follows 

O-LO = /lO ® fhO 
C^NLO = /nLO ® /nLO ® [1 + a5,NLO-f^*-^^] 
0-NNLO = /nNLO ® /nNLO ® [1 + a^.NNLO^^^^ + ("S.NNLo)^^^^^] (22) 



^For a precise measurement, we should allow for W~^W pair production and for W bosons produced via 
t-quark decays. These can contribute about 1% of the total signal. 
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where the label on as indicates the order to which the /3-function is evaluated. The parton 
distributions / obtained from the new global LO, NLO and NNLO analyses |]26[ correspond to 



as{Ml) = 0.1253, 0.1175 and 0.1161 respectively. The NLO and NNLO coefficient functions 
K^'^ are known |^. However, although the relevant deep inelastic coefficient functions are 
also known, there is only partial information on the NNLO splitting functions. Recently van 
Neerven and Vogt have constructed compact analytic expressions for the NNLO splitting 
functions which represent the fastest and the slowest possible evolution that is consistent with 
the existing partial information. These expressions have been used to perform global parton 
analyses at NNLO p6[. The uncertainty in the NNLO predictions of aw and az due to the 



residual ambiguity in the splitting functions is shown by the width of the NNLO bands in Fig. 4. 
This amounts to about ±1% uncertainty at the LHC energy, and less at the Tevatron. For 
completeness, we note that the dashed lines in Fig. 4 correspond to the quasi-NLO prediction 

c^NLO' = /nlo ® /nlo ® [1 + as,NLoK^^^ + (a5,NLo)^^^^^], (23) 

which was the best that could be done before the work of Refs. |2^, |26| . 

We see that the LO — > NLO — > NNLO convergence of the predictions for aw,z is good. The 
jump from ctlo to o"nlo is mainly due to the well-known, large, 0(0:5) vr^-enhanced Drell-Yan 
K-factor correction, arising from soft-gluon emission. The NLO and NNLO cross sections are 
much closer and, if this was the end of the story, W and Z production can clearly be predicted 
with sufficient accuracy. 

However let us turn to Figs. 5 and 6, each of which combine results presented in Refs. 



and [Q. The solid squares and triangles show the additional uncertainty in the predictions 
for aw and az which arise from changing the input information in the global parton analyses. 
The two major uncertainties appear to be due to the value of as and to using different parton 
densities labelled by qj and q|. The plots show the change in aw and az which is caused 
by changing the value of a5(M|) by ±0.005 respectively. The change in aw and az at the 
LHC energy is enhanced as compared to that at the Tevatron, since DGLAP evolution is more 
rapid at the smaller x values, x ~ Mw,z/\/s, probed at the higher energy. However, with our 
present knowledge of 0:5, the uncertainty Aas = ±0.005 is too conservative, and ±2% is a 
more realistic uncertainty in aw,z from this source at the LHC energy. 

The normalisation of the input data used in the global parton analyses is another source 
of uncertainty in aw,z- The HERA experiments provide almost all of the data used in the 
global analyses in the relevant small x domain. The quoted normalisation uncertainties of the 
measurements of the proton structure function F2 from the HI and ZEUS experiments vary 
with Q^, but a mean value of ±2.5% is appropriate. The qj and q| parton sets correspond to 
separate global fits in which the HERA data have been renormalised by ±2.5% respectively. 
For W and Z production, of qq origin, we naively would expect this to translate into a ±5% 
variation in aw,z, but the effect of DGLAP evolution up to ~ Mw,z is to suppress the 
difference in the predictions. 
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In summary, allowing for all the above uncertainties, we conclude that the cross sections of 
W and Z production are known to ±4% at the LHC energy, and to ±3% at the Tevatron. A 
major contributer to this error is the uncertainty in the overall normalisation of the HI and 
ZEUS measurements of F2. The normalisation may be made more precise by experiments in 
Run II at the Tevatron. 



5 Parton-parton luminosity 

In some circumstances it is sufficient to know the parton-parton luminosity, and not the proton- 



proton luminosity, see, for example, |T^. Of course if the proton-proton luminosity is known 
then the parton-parton luminosities can be calculated from the parton distributions determined 
in the global parton analyses. However in this case we rely on the normalisation of experiments 
at previous accelerators which yielded data that were used in the global analyses. 

Thus it may be better to monitor the parton-parton luminosities directly in terms of a 
subprocess which can be predicted theoretically to high precision. The best example is inclusive 



W (or Z) boson production, which is predicted up to two-loops, that is to NNLO The 
accurate observation of W (or Z) production at the LHC may therefore be used to determine 
the quark- ant iquark luminosity^]. Then the gluon flux, for example, may be determined from 
the global parton analysis which already is made to describe the measured W (or Z) cross 
sections. One advantage of this technique is that, for most LHC applications such as the search 
for new heavy particles, we need to consider a smaller interval of DGLAP evolution than has 
been the practice hitherto. 

Other ways to constrain the gluon-gluon luminosity are to study tt production or the pro- 
duction of two large pt prompt photons. The leading order subprocess is qq — > 77. However at 
LHC energies the gg — >■ 77 box diagram gives an important contribution. The two diagrams are 
shown in Fig. 7. The relative contributions are shown in Fig. 8 as a function of the transverse 
momentum of the photon pair = q^^ + q^2t the case when each photon has transverse 
momentum q^t > 20 GeV and rapidity jr^^l < 1. We see that in this kinematic domain the 
99 ~^ 77 subprocess gives a major contribution due to the higher gg luminosity. However 
there is a strong possibility of contamination by the subprocess gq 'j'jq, unless severe photon 
isolation cuts are imposed, see, for example, Ref. 

So far we have considered conventional parton distributions a{x, /i^) integrated over the 
parton transverse momentum kf up to the factorization scale fi. However many reactions are 
described by unintegrated distributions fa{x,k^,fi'^) which depend on both kt and the longi- 
tudinal momentum fraction x carried by the parton. In principle, unintegrated distributions 

""At first sight the qq luminosity in a given xi, X2 bin may be obtained by observing the number of W events 
in that bin and dividing by the qq ^ W cross section. However at NLO we include qq — > Wg^ qg — s- Wq etc., 
so the only possibility is to use the new W data to determine the qq luminosity within a global parton analysis. 
To do the same at NNLO we would require the NNLO expression for da/dy^, where yi is the rapidity of the 
decay lepton. 
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are necessary for the description of all processes which are not totally inclusive. In these cases, 
instead of the conventional QCD factorization, we have kt factorization Indeed, for some 
processes, after the specific integration over kt, the conventional 'hard' QCD factorization may 
be destroyed. 

The unintegrated quark-quark luminosity may be determined, for example, by observing 
the kt distribution of W (or Z) bosons or Drell-Yan pairs. At leading order, the fcrdependence 
of cross sections is given by the convolution 



da 
IPkt 



oc 



pi 



pI 



+P2 - kt) /g(a;i,p?,/i^) fq{x2,pl, fi^). 



(24) 



At large kt ^ Aqcd? to leading In kt order, the dominant contribution comes from either the 
domain pi ~ fct, P2 ^ h or the domain p2 ~ h, pi -C kt. It is natural to choose the 
factorization scale fi ^ kt. Then the kt distribution takes the simple form ||3^ 

da 



d^kt 

since at leading order 



oc 



xiqi{xi,kt) fg{x2,kt,kt) + X2q2ix2, fg{xi, kt , kt) 



dp^ 

o Jl 



a\X, p , kt) 



(25) 



(26) 



Moreover diffractive processes (with one or more rapidity gaps) are described, in general, by 
skewed parton distributions. An example is the double-diffractive exclusive Higgs boson pro- 
duction, described by the Feynman diagram shown in Fig. 9. This process is described by 
skewed gluon distributions with xi ^ x\ and X2 7^ x'2 [§3], ^ . 

Sometimes even the "skewed" parton flux can be monitored directly via another process with 
similar kinematics. For example, we may monitor the effective Pomeron-Pomeron luminosity 
for Higgs production of Fig. 9 by measuring double- diffractive dijet production in the region in 
which the transverse energy of the jets [Et) is of about half that of the Higgs mass [Et ~ Mn/'i) 

mm- 

Fortunately most of the physically relevant processes are described by skewed distributions 
with -C 1. In these cases the skewed distributions may be reliably reconstructed from the 
known conventional parton distributions [HO]. 



6 Conclusions 

We have studied the theoretical accuracy of the main three proposals for determining the 
luminosity of the LHC pp collider — namely using forward elastic data, lepton-pair production 
and W or Z boson production. The desired goal of a measurement to within ±5% seems 
theoretically attainable. 

We focused on potential shortcomings of each method. First, we demonstrated that the 
t dependence of the elastic cross section is well under control and, in fact, it turns out that 



12 



it can be safely approximated by a simple exponential in the region \t\ < 0.05 GeV^. For 

lepton-pair production we evaluated the corrections to the cross section arising from the strong 
interactions between the protons. We showed that in the relevant kinematic domain, with small 
transverse momentum of the produced lepton-pair, these effects are negligible. Hence a pure 
QED calculation of the cross section will give sufficient accuracy. 

The cross sections for W and Z production can now be predicted to NNLO, which at 
first sight would seem to provide an LHC luminometer with ±1% accuracy, see Fig. 4. How- 
ever the uncertainties in the input to the global parton analyses mean that the error could, 
conservatively, be as large as ±4%. The uncertainty may be reduced if we work in terms of 
quark- antiquark luminosity, which is relevant in some future applications of the LHC. 

Finally, we note that luminosity determinations based on the measurement of the forward 
elastic cross section and (most probably) on two-photon e+e~ production can only be made 
in low luminosity runs, and require dedicated detectors and triggers. On the other hand, the 
measurement oiW or Z and two-photon production may be performed at high luminosity 
with the central detector and with standard triggers. 
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(a) 



(b) 



Figure 1: (a) production mediated by photon exchange, and (b) possible rescattering 
corrections. 
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Figure 2: (a) production mediated by photon exchanges and (b) two of the possible 
rescattering corrections. 
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Figure 3: (a) Lepton pair production in pp collisions, (b) one of the rescattering corrections, 
and (c) a possible contamination coming from proton dissociation into X, Y systems. 
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Figure 4: The predictions of the cross sections for W and Z production and leptonic decay at 
the Tevatron and the LHC obtained from parton sets of LO, NLO and NNLO global analyses 
of the same data set |2^. The cross sections labelled LO, NLO and NNLO are as in (|2^ ) and 
the dashed line is the NLO' prediction of (|2^) . The band of the NNLO predictions allow for 
the ambiguity in the NNLO splitting functions p8|. Also shown are measurements obtained at 
the Tevatron 



30[| . The figure is taken from p6 |. 
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W and Z Cross Sections: Tevatron 
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Figure 5: The solid squares and triangles are the predictions of the NLO' cross sections of 
( |25| ) for W and Z production and leptonic decay in pp coUisions at a/s = 1.8 TeV obtained 
using various NLO sets of MRST99 partons The open square and small error bar are, 

respectively, the NLO and NNLO predictions of (|22D using the MRSTOO partons [||]. Also 
shown are the experimental measurements from CDF and DO For ease of reference 
±5% lines are shown about the MRST99 default prediction. We thank W.J. Stirling for this 
figure, which combines results presented in pSl p6|. 
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W and Z Cross Sections: LHC 
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Figure 6: As for Fig. 5 but for pp collisions at = 14 TeV. 
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Figure 7: Feynman diagrams driving 77 production in pp collisions. 
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Figure 8: Contributions to the diphoton pt spectrum from the diagrams of Fig. 7 in pp coUisions 
at ^/s — 14 TeV. Each photon is required to have transverse momentum pjt > 20 GeV and 
rapidity \r]^\ < 1. The photons are required not to he within the same rj — (p cone of radius 0.4. 
We thank M.A. Kimber for this figure. 
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Figure 9: Higgs production via Pomeron-Pomeron fusion in pp collisions. 
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